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AN IMPROVED ANALYTICAL TREATMENT OF THE DENTING 
OF THIN SHEETS BY HAIL 


By Robert G. Thomson and Robert J. Hayduk 
Langley Research Center 

SUMMARY 

Severe structural damage, such as surface erosion, dents, perforations, and tears, 
can occur when high-speed aircraft collide with hailstones. In this paper the denting of 
aircraft skin by hail is mathematically analyzed and the results are compared with 
experimental data. The denting process is modeled mathematically by assuming that a 
crushable, spherical hailstone impacts normal to a flat sheet. The resultant dent depth 
and shape are determined by utilizing an existing computer program which considers 
both bending and membrane action and elastic and plastic material behavior. The results 
of this analysis are compared with experimental data from the British Royal Aircraft 
Establishment and with a previous analysis which considered bending only. The improved 
analytical treatment of the denting process agrees well with the experimental data and 
shows that membrane forces must be considered when sheet deflections are large. 

INTRODUCTION 

Severe structural damage can occur when high-speed aircraft collide with hail- 
stones. Large hailstones erode surfaces, cause dents, and perforate and tear sheet- 
metal skins. References 1, 2, and 3 show photographic evidence of such damage to 
leading edges, empennages, and forward sections as well as photographs of shattered 
windshields. Methods of predicting hail damage to aircraft surfaces are needed to 
determine how aircraft should be designed to survive an accidental encounter with hail. 

The skin deformations resulting from hail impact usually are plastic near the point 
of impact and elastic in the surrounding area. Reference 4 was a first attempt at pre- 
dicting the denting characteristics of hail on thin sheets. In that reference the permanent 
deformations were predicted by an approximate elastic-plastic analysis. Classical thin- 
plate bending theory was used to circumvent the difficulties of a large-deformation anal- 
ysis which would include both bending and membrane effects. However, recent hail- 
impact data from the British Royal Aircraft Establishment (R.A.E.) indicate that 
membrane forces must be included. (See ref. 5.) The results of the bending analysis of 
reference 4 deviate from these experimental data when sheet deflections are comparable 


with sheet thickness or larger. A similar conclusion was reached by A. L. Florence in 
reference 6, which presents the results of an experimental investigation of circular plates 
subjected to uniformly distributed impulses. Florence compared the experimentally 
obtained permanent deformations with those predicted by rigid-plastic bending theory 
and concluded that a "full treatment of the problem taking account of the membrane 
action is certainly required.” 

The evidence indicates that large-deformation theory is required to analyze the 
hail-impact problem; therefore, this paper presents an improved analytical treatment 
which includes the effects of membrane action. The solution was accomplished by 
adapting the DEPROSS 3 computer program to the hail-impact problem. DEPROSS 3 
was developed at Massachusetts Institute of Technology by E. A. Witmer, H. A. Balmer, 

J. W. Leech, T. H. H. Plan, and W. Herrmann (see refs. 7, 8, and 9) to analyze the 
dynamic elastic-plastic response of impulsively loaded simple structures. With minor 
modifications in the initial velocity input section and the addition of the hailstone mass to 
the mass of the plate, the program was employed to predict the size of the dents caused 
by hail. 

A discussion of the modifications made to DEPROSS 3, the final permanent defor- 
mations obtained, and comparisons with the bending solution of reference 4 and the R.A.E. 
experimental data are presented herein. 

The authors gratefully acknowledge the assistance offered them by the British 
Royal Aircraft Establishment. Through personal communication with Ian I. McNaughton 
of the R.A.E. , the authors received supplementary data on hail impact with aircraft. 

These supplementary data included unpublished results for 1.27- and 1.91-cm- (0.5- and 
0.75-inch-) diameter hail-impact for flat sheets. 

SYMBOLS 

The units used for the physical quantities defined in this paper are given both in the 
International System of Units (SI) and in the U.S. Customary Units. (See ref. 10.) Appen- 
dix A presents factors relating these two systems of units. The measurements and cal- 
culations were made in the U.S. Customary Units. 

A radius of crushed hailstone, aK 2 

a radius of hailstone 

b radius of central plastic region, where Mj. = 0 
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c 



D flexural rigidity of sheet, Eh^/l2(l - t>2) 

E elastic modulus of sheet material in tension and compression 

gp hailstone impact velocity 

h thickness of sheet 

i integer 

Jn(z) Bessel function of the first kind, where n is the order of the Bessel function 

Kl factor which adjusts the height of the crushed hailstone, 



K 2 factor representing the instantaneous spread of the base of the hailstone 

(i.e., a factor of 1.5 or 2.0) 

k = - 2^2) 

Mq yield-moment resultant, aQ(h/2) ; maximum bending moment that sheet can 

sustain 

Mr,Mg radial and circumferential bending-moment resultants 
m mass distribution 

p Hankel transform parameter 

R radius of finite sheet 

r,z radial and transverse coordinates (see fig. 1) 

ASi distance between mass stations, where i = 1, 2, . . . 
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t 


time 


t* time of cessation of all sheet motion after impact 

KiA Pp 
— 

V = ° 

° ^ ^ KiA Pp 

h Pg 


(ki = 


2 for bending-only analysis (see ref, 


4 )) 


* 


w sheet deflection in transverse direction (z -direction) 

6 permanent center deflection, w(0,t*) 

rj - r/A 

jjL mass per unit area of sheet, Pghg 

V Poisson's ratio 



p mass density 

Oq yield stress in simple tension or compression (assumed to be identical in 

magnitude) 

T = t/k 

momentum constant, set equal to 0.70 (see ref. 4) 

Subscripts: 


p projectile 

s sheet 

A dot over a symbol indicates partial differentiation with respect to time t. 
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COMPUTER ANALYSIS 


Program Description 

The DEPROSS 3 computer program used in this analysis is described in detail in 
references 7, 8, and 9. A brief description of the program will be included in this sec- 
tion, but the emphasis will be on the minor modifications made to the program and the 
results obtained from it. The modifications adapted the program to the specific problem 
of hail impact on thin sheets. 

DEPROSS 3 is a numerical computer program capable of predicting the dynamic, 
axisymmetric response of shells, plates, rings, and beams to impulsive or blast loads. 
Material behavior can be elastic, elastic-perfectly-plastic, elastic-strain-hardening, or 
elastic-strain-hardening with strain-rate sensitivity. The elastic-strain-hardening 
material behavior was chosen to represent the response of the thin sheets to hail impact. 
A bilinear stress-strain relationship was adopted in which the target material was 
assumed to behave elastically up to the yield stress at 0.2-percent offset and then to 
experience linear strain-hardening between the yield stress and the ultimate stress. The 
coordinates of the bilinear stress-strain curve used in the computer program are listed 
in table I. 

The model which represents the simple structure consists of concentrated masses 
connected by massless links. The links can extend or contract axially but are infinitely 
rigid in bending; all bending action takes place at the mass points. External forces are 
applied at these mass points; internal forces and moments are transmitted between the 
masses by the links. The dynamic equilibrium equations are solved numerically in 
finite-difference form for each mass point at successive time increments. The incre- 
ments in the forces and moments at each mass location are determined by representing 
each mass station as consisting of a number of discrete layers of material. Each layer 
carries normal stresses and is separated from an adjacent layer by a material that has 
infinite shear rigidity and does not carry normal stresses. The axial force and bendii^ 
moment at each mass location can then be determined from the normal stresses in the 
various layers. The strain in each layer is determined by adding the midplane strain and 
the bending strain due to curvature, that is, plane sections remain plane during deforma- 
tion. Rotary inertia and transverse shear deformations are assumed to be negligible. 

Determination of the Mass Distribution for the Program 

The flat surface of a wing panel is represented in the DEPROSS 3 program as a flat, 
circular, clamped sheet of uniform thickness. The radius of the circular sheet is taken 
to be large enough so that the effects of the clamped boundary on the central deflection of 
the sheet are negligible. In reference 5, it was observed that "the indentation base shape 
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(of the hailston^ is formed very rapidly and is maintained throughout the impact phase 
being progressively driven deeper through the plate. In the case of 1 inch hail impacting 
on L73 material [2014-T6, U.S. equivalent , the base shape is that of a 2 inch diameter 
sphere and on L72 material [2014-T4, U.S. equivalent it is a 1.5 inch diameter sphere." 
Two constants and K 2 ^ were introduced into the analysis to account for this experi- 
mental observation. The constant K 2 in the present analysis represents the instanta- 
neous spread of the base of the hailstone (i.e., a factor of 1.5 or 2.0). The height of the 
crushed hailstone was adjusted by the factor Kj 


Ki = 



( 1 ) 


to yield the same mass as the original uncrushed hailstone. 

The mass distribution of this crushed hailstone and sheet is shown in figure 1(a). 
The radial coordinate r of the sheet is nondimensionalized by dividing by the base 
radius of the crushed hailstone, A = aK 2 , 


The sheet material was concentrated in 32 mass stations equally spaced along the radius, 
as shown in figure 1(b). This number of mass stations adequately represented the 
dynamic response of the sheet. (See ref. 8 .) The first and last mass stations were one- 
half of a space away from the origin and boundary, respectively. The mass of the crushed 
hailstone was also distributed among the mass stations. The mass stations were chosen 
so that an integer number of mass points were beneath the hailstone. The mass distribu- 
tion of the crushed hailstone is given by 


mp(T?) = KiA 



(OSrjSl) (3) 


This exact mass distribution is shown in figure 1(b) along with the stepwise approxima- 
tion. The additional mass contributed by the hailstone at each station was calculated by 
taking the mass distribution as uniform between stations and equal to the value obtained 
from equation (3) for the 77 's corresponding to the mass stations. The mass at the last 
station adjacent to the periphery of the hailstone was adjusted so that the total mass 
included in the stepwise DEPROSS 3 distribution was equal to the total mass of the 
hailstone. 


6 



In the calculations, the number of mass stations beneath the projectile varied with 
hail size and target material. When the hail diameters were modified by the K 2 factor 
of 1.5 and 2.0 for L72 and L73, respectively, there were 6 or 8 mass stations beneath the 
2.54-cm- (1-inch-) diameter hail, 4 or 6 beneath the 1.9 1-cm- (0.75-inch-) diameter 
hail, and 3 or 4 beneath the 1.27-cm- (0.50-inch-) diameter hail. For the 1.91-cm- 
(0. 75-inch-) diameter hail and L72 target material, the plate radius was increased from 
10.2 to 11.4 cm (4 to 4.5 inches) so that an integer number of mass stations was com- 
pletely covered by the hailstone. An integer number of mass stations ensured a smooth 
velocity distribution. The details of the programing are given in appendix B, statement 
numbers 8 through 104. The computer program printed the displacements of the mass 
stations every 0.06 iis. 

Determination of the Initial Velocity Distribution for the Program 

The initial velocity distribution of the sheet and crushed hailstone was determined 
by assuming that the momentum exchange upon impact was instantaneous and axisym- 
metric. The momentum balance for a completely inelastic impact is 


2ttA2t 7 dT7[mp(77) + = mp(77)go277A2T7 dr? (4a) 

or 

2 irA ^?7 drj^iA^jl - V^'jpp + 2 ttA^t] dtj hp^w(? 7 , 0 ) = 27tA^T 7 d77|KjA^l - V^^Pp% 
Solution of equation (4b) for the initial velocity distribution gives 


w(tj,0) 


^1^ Pp 
h Ps 




KiA 

1 +-i- 

h 



(0 ^ T? ^ 1) (5a) 


w(?7,0) = 0 (77 = 1) (5b) 

The initial velocity distribution (eqs. (5)) for the combined mass of hailstone and 
sheet is illustrated in figure 1(c). The initial velocity imparted to each mass station was 
calculated from equations (5) using the 77 's corresponding to the stations except for the 
last mass station adjacent to the outer periphery of the hailstone. The velocity of this 
last mass station was adjusted so that the integrated, stepwise, initial velocity distribution 


L_ 
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(see fig. 1(b)) was the same as the integrated, exact velocity distribution, equations (5). 
The program listing for the initial velocity distribution is given in appendix B (statement 
numbers 50 through 85). 

EXPERIMENTAL DATA FROM BRITISH ROYAL AIRCRAFT ESTABLISHMENT 

The R.A.E. has conducted extensive hailstone impact tests in order to obtain basic 
data on the hail-impact resistance of aircraft. The results of an experimental program 
using flat aluminum sheets as targets and 2.54-cm- (1-inch-) diameter hail are pre- 
sented in reference 5. Additional data for other hailstone sizes were personally com- 
municated to the authors. These supplementary R.A.E. data were for 1.27-cm- 
(0. 50-inch-) diameter and 1.91-cm- (0.75-inch-) diameter hail. The target sheets were 
made from two different British aluminum alloys, designated L72 (2014-T4, U.S. equiv- 
alent) and L73 (2014-T6, U.S. equivalent). The mechanical properties of L72 and L73 are 
listed in table I. The sheets were 30.5 cm (12 inches) square and were bolted to a frame 
leaving a 20-cm- (8-inch-) square, unsupported area in the center. The sheet thick- 
nesses were 0.071 cm (0.028 inch), 0.091 cm (0.036 inch), 0.122 cm (0.048 inch), and 
0.163 cm (0.064 inch). 


RESULTS AND DISCUSSION 

The modified DEPROSS 3 program was used to calculate the permanent deforma- 
tions of thin, flat sheets impacted by hailstones of various sizes. The sheet thicknesses 
and hailstone sizes corresponded to thicknesses and sizes used by the R.A.E. in its 
experimental studies. Since DEPROSS 3 handles only axisymmetric problems, the 
square experimental targets were idealized as 10.2-cm- (4-inch-) radius clamped cir- 
cular plates. The difference in the shape between the experimental and idealized targets 
seems to be unimportant inasmuch as the permanent depth of the dent was established by 
calculations before the disturbance reached the boundary of the target. Additional cal- 
culations were made using the classical bending approach of reference 4 for comparison 
purposes. The following sections will discuss these calculations and comparisons. 

Sheet Displacements 

The computed shapes of the sheet at various times (24, 102, 204, and 420 fxs) are 
shown to scale in figure 2 for two impact conditions. After impact the disturbance travels 
radially out to the boundaries as the plate deflects. The deflection of the central region 
becomes quite stable before the disturbance reaches the clamped boundaries. Since 
DEPROSS 3 does not include internal damping, the plate continues to oscillate about its 
plastically deformed position. 
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In figure 3 the center deflection (taken to be the displacement of mass station 1) of 
a 0.091-cm- (0.036-inch-) thick L72 plate is plotted as a function of time for 1.27-cm- 
(0. 50-inch-) diameter hail, which impacted at various velocities. The center deflection 
increases with time to a maximum and then oscillates elastically. Additional calculations 
indicate that oscillations are more pronounced for the larger hail sizes. After approxi- 
mately 300 jus, plastic work ceases and the sheet has absorbed approximately 82 percent 
of the initial kinetic energy of the hailstone. Twelve percent of the remaining energy is 
elastic -strain energy, and the other 6 percent is kinetic energy. A mean value of the 
center deflection after 300 /is is taken as the permanent center deflection. 

After long periods of computation, instabilities in some of the numerical calcula- 
tions occurred. These instabilities are not easily recognized until they become pro- 
nounced. In figure 3, the computations for g^ equal to 66 m/s (216 ft/s) and 183 m/s 
(600 ft/s) are suspected of developing numerical instability. These instabilities were 
minimized by carefully choosing the time increments between calculation cycles. (See 
ref. 9.) 

Figure 4 shows comparisons between the computer-predicted final shapes and the 
R.A.E. experimentally determined shapes for L72 sheets. The hailstones were 1.27 cm 
(0.50 inch), 2.54 cm (1 inch), and 1.91 cm (0.75 inch) in diameter (represented in the 
figures by their radii). There is a close correlation between the experimental data and 
the computer predictions; this correlation is much closer than the conical shapes pre- 
dicted by the classical bending approach of reference 4. The computer predictions 
exhibit slightly larger overall deflections than do the experimental data and, hence, are 
conservative. This overprediction can be partially attributed to the neglect of strain- 
rate effects (see ref. 11, p. 64) in the calculations. Previous comparisons between 
DEPROSS 3 calculations and experimental data also indicated a slight conservatism in 
the calculations. (See ref. 8.) Reference 8 indicated that the specific factors responsible 
for this overprediction had not been isolated and resolved, but among the possible factors 
were an ". . . inadequate representation of the stress-strain properties of the material, 
neglect of strain-rate effects, an inadequate representation of the flow process, uncer- 
tainties of the externally applied forces, etc." 

The lateral splash of hail on impact, observed experimentally by the R.A.E. , could 
also contribute to this conservatism. After impact, an annulus of material from the 
crushed hailstone continues to move radially outward. This radial dispersion of mass 
is not accounted for in the DEPROSS 3 calculations. Another possible factor is the dif- 
ference between the experimental sheets and computer idealization of those sheets. The 
experimental sheets were 20.32 cm (8 inches) square, whereas a 20.32-cm- (8-inch-) 
diameter circular plate was used in DEPROSS 3. The square sheet contained 27 percent 
more material area than the circular plate; however, this additional area was near the 
boundary and should have little effect on the permanent center deflection. 
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Permanent Center Deflection 


In figures 5, 6, and 7, the computer predictions are compared with the experimen- 
tally determined permanent center deflections of L72 aluminum. The center deflections 
for three different sheet thicknesses are shown in each figure. 

The supplementary R.A.E. data for L72 sheets impacted by 1.27-cm- (0.50-inch-) 
diameter hail are presented in figure 5. A solid line has been faired through the data, 
and the dashed line represents the theory. The theory generally overestimates the 
deflection. A vertical bar marked with either P or SP indicates, respectively, that 
the sheet was either completely penetrated or developed a split at the corresponding 
velocity. Supplementary data for 1.91-cm- (0.75-inch-) diameter hail are shown in fig- 
ure 6. Again, the theory is generally conservative, but still closely follows the slope of 
the experimental data. 

Figure 7 contains the data for 2.54-cm- (1-inch-) diameter hail taken from refer- 
ence 5. Again, the theory is conservative but follows the slope of the experimental data. 

Also shown in figure 7, by the three solid symbols, are computer results for a 
simply supported circular plate with the same dimensions and loading conditions as the 
corresponding clamped circular plates. Calculations were made for two simply supported 
plates: the0.122-cm- (0.048-inch-) thick plate and the 0.091-cm- (0.036-inch-) thick 
plate. These results show that the boundary conditions of the plate have little effect on 
the permanent center deflection. These calculations for the simply supported plates also 
exhibited a lower frequency of elastic oscillation, as might be expected. Consequently, 
longer computation times were needed to establish the permanent deformations. 

Comparison of Bending and Bending-Membrane Theories With 

Experimental Data 

A comparison of the classical, thin-plate bending analysis of reference 4 and the 
large -deflection bending and membrane numerical analysis (DEPROSS 3) with the R.A.E. 
experimental data is shown in figures 8 and 9 for L72 and L73, respectively. In the fig- 
ures the permanent center deflection of a 0.122-cm- (0.048-inch-) thick sheet is plotted 
as a function of hailstone impact velocity. Curves are presented in each figure for three 
different hailstone sizes, 1.27-, 1.91-, and 2.54-cm (0.50-, 0.75-, and 1-inch) diameter 
(represented in the figures by their radii). 

The bending analysis of reference 4 was modified to include the experimentally 
observed spreading of the hailstone. This modification was accomplished by increasing 
the hailstone diameter by the factor K2 discussed previously and decreasing the hail- 
stone density to maintain the original mass. It was not possible to include the mass of 
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the hailstone after impact in the bending analysis without resorting to major modifica- 
tions. The appropriate equations and the details of the calculations are contained in 
appendix C. 

The bending-membrane theory follows the trend of the experimental data quite well, 
while the bending theory deviates from the data as soon as the deflection increases above 
the sheet thickness. The experimental data indicate a linear dependence on the hailstone 
velocity g„ whereas the bending theory of reference 4 yields too strong a velocity 

o Oc 

dependence, that is, 6 cc g^- . The bending-membrane results for L72 follow the linear 

slope of the experimental data but are conservative. The results from the bending- 
membrane theory for L73 show very good agreement with the experimental data both in 
slope and magnitude. The better agreement with L73 data than with L72 data is probably 
due to a better approximation of the stress-strain behavior of L73. There may also be a 
strain-rate effect, but data on the effects of high strain rates on the behavior of these 
materials were not available. 


Empirical Relationship 

An attempt was made to correlate the computed center deflection curves of fig- 
ures 5, 6, and 7 with one parameter. Logarithmic cross plots of the calculations were 
made to determine exponents, and the result of these efforts is shown in figure 10 where 
the center deflection is plotted as a function of gQa2y4i3/4 on log-log scale. The 
parameter reasonably correlates the curves. It is apparent that the exponent on "a" is 
not a constant since the curves shift as the value of "a” changes. A functional relation- 
ship more complex than "a" raised to a power is required to better correlate the curves. 


CONCLUDING REMARKS 


In this paper the denting process which occurs when a high-speed aircraft collides 
with hailstones was modeled mathematically by assuming that a crushable, spherical 
hailstone impacted normal to a flat sheet. The solution was accomplished by adapting 
the DEPROSS 3 computer program to the hail-impact problem. This analytical method 
considers both bending and membrane action and elastic and plastic material behavior. 

The analytical results and experimental data from the British Royal Aircraft 
Establishment presented in this paper suggest that membrane action is a dominant fac- 
tor in the denting of aircraft surfaces by hail when the deflection exceeds the sheet thick- 
ness. According to reference 4, the center deflection should vary with the hailstone 
velocity to the 2.85 power when membrane forces are neglected. The improved analysis 
of this paper, which includes membrane action, indicates that the center deflection varies 
with velocity to the first power. 
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The bending and membrane analysis of DEPROSS 3 agrees well with the experi- 
mental data for various combinations of hailstone sizes and aluminum -sheet thicknesses. 
Although approximate, the bilinear stress -strain relationship employed in the computer 
program gave results that quite closely paralleled the trend of the experimental data. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va. , December 4, 1970. 
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APPENDIX A 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (SI) was adopted by the Eleventh General 
Conference on Weights and Measures, Paris, October 1960, in Resolution No. 12 (ref. 10). 
Conversion factors for the units used herein are given in the following table: 


Physical quantity 

U.S. Customary 
Unit 

Conversion 

factor 

(a) 

SI Unit 

Length 

in. 

0.0254 

meters (m) 

Density 

Ib-s2/ft4 

16.02 

kilograms per cubic meter (kg/m^) 

Stress 

psi (lbf/in2) 

6.895 X 10^ 

newtons per square meter (N/m2) 

Velocity 

ft/s 

0.3048 

meters per second (m/s) 


^Multiply value given in U.S. Customary Unit by conversion factor to obtain 
equivalent value in SI Unit. 


Prefixes to indicate multiple of units are as follows: 


Prefix 

Multiple 

micro (p) 

10-6 

centi (c) 

10-2 

kilo (k) 

102 

giga (G) 

109 
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APPENDIX B 


MODIFICATIONS TO DEPROSS 3 

A listing of the velocity input section to the modified DEPROSS 3 computer program 
is presented in this appendix. The four subroutines PRINT, EQIHL, STRAIN, and STRESS 
remain essentially the same as presented in reference 9. The lack of identification num- 
bers at the right side of the following listing indicates modifications made to the original 
DEPROSS 3 velocity input section. 

Main Program 

PROGRAM DEPROS3 ( 1 NPuT . OUTPUT « TAPE5= I NPUT « TAPES ) 

C CLAMPED CIRCULAR PLATE. AXISYMMETRIC MOTION (DEPROSS 3) 

C INITIAL VELOCITY DISTRIBUTION 

DIMENSION R (50 ) .DR (50 ) . Z (50 ) .DZ (SO 1 . DELTS (50 ) . DOS ( 50 ) .DELTT (50 ) . DD 
I T (50 ) .DT (50 ) .S INT (50 ) . DSINT (50) .COST (50 ) .EPS I TI (50 ) .EPS I TO (50 ) .EPS 
2 IP I ( 50 ) .EPS I PO ( 50 ) .8 I GNT ( 50 ) .BIGNTR (50 ) .8 I GNP ( 50 ) . QTR (50 ) .B I GMT ( 50 
3 ) .BIGMP(50 ) .RO (50 ) .ROAVE (50) .CRO (50 ) . S I NTO ( 50 ) . CMASS ( 50 ) . Y I ELD ( 50 ) 

4 .ZETA (5 ) .CZETA (5 ) .RZpTA (5.50 ) .E (6 ) .EPSIL (5 ) . S IGMA ( 5 ) . SNOSQ (5 ) .HSFL 
5 (5 > . SNT (5. 5. 50 ) .SNP ( 5 . 5. 50). S( 50) , VELOCY (50) .N5TAT(50) .TH(50) .C5(5 


60 ) 

COMMON R. DR. Z.DZ. DELTS .DDS .DELTT .DDT .DT , S INT . DS INT . COST , EPS I T I . EPS 000 7 
I I TO.EPS IP I . EPS IPO .81 gNT . B I GNTR . B I GNP . QTR . B IGMT .3 1 GMP.R0 .ROAVE . CRO. 0008 
2S INTO. CMASS. Y I ELD. ZETA .CZETA .RZETA . E.EPS IL . SIGMA . SNOSQ. HSFL . SNT . SN 0009 
3P.N.NI .N2 .NFL .NSFL .NYl ELD.MY lELD.KE . J . T I ME .DELTSO . SNU , C I NETO . HHALF 0010 


4.C2.C3.C4.C6.C7.M 
5000 CONTINUE 

READ I 0 . N. NFL .NSFL . M , M 1 . M2 . KE .8 I GR . H . CFL . SNU . RHO . DEL TA T . D . P . (EPSIL 
1 (L > .SIGMA (L ) .L = l .NS^L ) 

IF (EOF .5 >400.40 1 

400 STOP 

401 CONTINUE 

READ 1 01 .8 IGA .RHP.GO. Cl 0 

C MV=0 A RESTART. MV=NEGATIVE NUMBER-STOP. Mv=POSITIVE NUMBER 
C -CALCULATES WITH INITIAL INPUT 
read 10. MV 
REWIND 8 

101 format (4E1 8.8 ) 

PRINT 9 .N.NFL .NSFL.M ,M 1 .M2 .KE . B I GR . H . CFL . SNU. RHO . DELTAT . D .P . ( EPS I L 
1 (L ) .SIGMA (L ) .L=l .NSFL ) 

9 FORMAT (//6X1 HN . 4X3HNFL . 3X4HNSFL . 5X 1 HM . 6X2HM 1 . 5X2HM2 . 5X2HKE/7 I 7/X9X 
14HB IGR. 1 5X1 HH. 1 6X3HCFL. I 5X3HSNUZ4E 1 8 . 8//9X3HRHO . I 4X6HOELTAT. 1 4X1 HD 
2. 1 7X1HP/4E18.8//8X5HFPSIL. I 3X5HS 1 GMA/ ( 2E 1 8 . 8 ) ) 

PRINT I 02. B IGA .RHP.GO. Cl 0 

1 02 FORMAT ( /Z9X4HB IGA . 1 4X3HRHP . 1 6X2HGO. 1 6X3HC 1 0/ (4E1 8.8 ) ) 


! 0 FORMA T(7I5Z(4E18.8) ) 

PRINT 8 . MV 
8 F0RMAT(//4H MV=15) 

Nl= N+1 0020 

N2= N 0021 

NYIELD= 0 (5022 

DELTS0= BIGRZFLOAT(N) 0023 

DELT= DELTS0Z2. 0024 
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APPENDIX B — Continued 


FLONFL= NFL 0035 
HFL= H/FLONFL 0026 
HHALF= H/2. 0027 
DHALF= HHALF/SQRT (FLONFL**2-CFL) 0028 
E(l)= SIGMA ( 1 )/EPSIL ( 1 ) 0029 
Cl= E( 1 )/{ 1 .-SNU**2 ) 0030 
C2= Cl/DELTSO 0031 
C3= -2.*SNU 0032 
C4* 3. 1A1592653590*DFLTS0/E( i » 0033 


DELTAR= BIGR/FLOAT (N ) 

Rl= DELTAR/2. 

LOAD= BIGA/DELTAR 
MOAD= LOAD-1 
IF (BIGA.GT.BIGR ) MOAO=N 
STH= 0. 

t C9=2./( (C10**3- (C10**2-l • )**1 .5)*C10**2> 

IF (MOAD.LE.O ) GO TO 3 

DO 3 NV=l*MOAD 

S<NV)= (2,*FL0AT{NV)-1 . )*R1 

TH (NV )= C9»SQRT (B1GA**2-S (NV )#»2 ) 

STH= (2.*FLOAT<NV)-l , )*TH(NV)+STH 
NSTAT<NV)= NV 

3 CONTINUE 

IF (BIGA.GT.BIGR ) GO TO 4 

NV» LOAD 

IF(NV.EO.O)NV=1 

S(NV)= (2.*FL0AT(NV)-1 . )*R1 

TH (NV >= B IGA **3*1 , /3 , *C9/ ( 2 . *R 1 **2* ( 2 . *FLOAT (NV )- 1 . ) l-STH/ (2.*FL0A 
1 T(NV )-l . ) 

NSTAT(NV)= NV 

4 CONTINUE 
Ll= LOAD 

IF(BIGA.GT.BIGR) Ll= N 

PRINT 104. (NSTAT(NV) , TH ( NV ) . NV= 1 « L 1 > 

1 04 FORMAT ( //3X5HNST AT . 9X2HTH/ ( I6.E18.6 )) 

DO 5 1=1 .LOAD 

C5 ( I )= DELTAT**2X (RH0*H*DELT50+RHP*TH ( I ) *DELTSO ) 

5 CONTINUE 
NOAD= LOAD+1 
DO 6 I=N0AD.N1 


C5(I)= OELTAT**2/ (RHO*H*DELTSO ) 

6 CONTINUE 

IF(D)210.12. 1 1 0035 

1 1 C6= 1 . /P 0036 

C7= 1 . 1547005/D/DELTAT/Cl 0037 

GO TO 13 0038 

12 C6= 1. 0039 

C7= 0. 0040 

13 IF (NSFL-1 >21 0 . 1 6. 1 4 0041 

14 DO 15 L«2.NSFL 0042 

E(L)= (SIGMA (L )-SIGMa (L-1 ))/ (EPSIL (L )-EPSIL (L-1 1 ) 0043 

15 CONTINUE 0044 

16 E(NSFL+n= 0. 0045 

DO 17 L=1 .NSFL 0046 

HSFL(L)= HFL* (E (L )-E (L+1 ) )/E ( 1 1 0047 

SN0SQ(i_)= (Ed )*EPSIl(L) )**2 0048 

17 CONTINUE 0049 

NFL1= NFL/2 C050 

DO 20 K=1 .NFLl 0051 


‘t’Cg = Kj (see eq. (1)). 
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APPENDIX B — Continued 


L= NFL-K+l 0052 

ZETA(K)» DHALF«FLOAT (L-K ) 0053 

ZFTA<L)= -ZETA(K) 0054 

CZFTA(K)» C2*2ETA(L) 0055 

20 CZETA<L)= -CZETA(K) 0056 

DO 21 1=1 .N1 0057 

SINTOC 11=0, 0058 

R0(I)= DELT+FLOAT n + I-1 ) 0059 

CR0(1)= C1/R0(I) 0060 

CMASS(I)= C5(1)/R0CI) 

21 CONTINUE 0062 

DO 23 I=2.N1 0063 

ROAVEt n = (R0( M+R0( I-l ) )/2. 0064 

Cl R=-C1/R0AVE ( I ) 0065 

DO 22 K=1 ,NFL 0066 

RZETA (K« I )=C1R*ZETA (K ) 0067 

22 CONTINUE 0068 

23 CONTINUE 0069 

IF(MV)210,90t30 0070 

30 MY1ELD= 0 0071 

J= O 0073 

TIME= 0. 0073 

DO 50 1*1 «N1 0074 

R( I )= R0( I ) 0075 

DR( I )= 0, 0076 

Z( I 1= 0. 0077 

DZ(I)= 0, 0078 

DELTS(1)= DELTSO 0079 

DELTTI I »= 0, 0080 

SINT( I )= 0. 0081 

COSTin* 1. 0082 

BIGNTI I )= 0. 0083 

BIGNPI I )= 0. 0084 

BIGMTt I )= 0. 0085 

BTGMP( 1 >* 0. 0086 

YIELD! I )= 555555550000000000008 0087 

DO 40 <=1 ,NFL 0088 

DO 39 L=1.NSFL 0089 

SNT(L»K« I )= 0. 0090 

SNP(L.K*n= 0. 0091 

39 CONTINUE 0092 

40 CONTINUE 0093 

50 CONTINUE 0094 


Cl= C9*BIGA/H*RHP/RH0 
DELTAR= BIGR/FLOAT(N) 

Rl= DELTAR/2. 

LOAD* BIGA/DELTAR 

noad= load-1 

IF(BIGA.GT«BIGR) MOAn=N 
SZDOTO* 0. 

IF (MOAD.LE.O ) GO TO 80 

DO 80 NV=l«MOAD 

S(NV)= (2,»FLOAT(NV)-l . )4R1 

VELOCY(NV)= C1*G0*SQRT{BIGA**2-S(NV)**2)/(BIGA+C1*SQRT(BIGA**2-S(N 
1 V )**2 ) ) 

SZDOTO= f2.*FLOAT (NV )-l . ) »VELOC Y ( NV ) +SZDOT0 
ZDOTO* VELOCY(NV) 

I 1 = NV 
12= NV 

NSTAT(NV)= NV 
DELTAZ= DELTAT*ZDOTO 
DO 70 1=11,12 
DZ ( 1 )= OELTAZ 
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APPENDIX B - Continued 


70 CONTINUE 0103 

80 CONTINUE 0104 

IF(BIGA.GT.BIGR) GO TO 86 
NV* LOAD 
IF<NV»EQ»0)NV«1 
S(NV)“ {a.*FLOAT(NV)-l . )*R1 

VELOCYINV )= BIGA**2*gO* < .5-1 ./Cl +1 ./Cl**2*ALOG ( 1 .+C1 ) )/ ( 2 . *R1 **2 * ( 
12.*FLOAT(NV)-l . ) )-(SZDOTO/(2.*FLOAT(NV)-l . ) ) 

ZDOTO* VELOCY(NV) 

11= NV 
12= NV 

NSTAT(NV)= NV 
DELTAZ= DELTAT*ZD0T0 
DO 85 1=11*12 
DZU)= DELTAZ 

85 CONTINUE 

86 OZ(N»= 0. 0105 

Ll“ LOAD 

IFfBIGA.GT.BIGR) Ll= N 

PRINT 103. (NSTAT (NV ) ,S (NV) « VEL0CY(NV ) *NV=1 *H ) 

1 03 FORMAT ( //3X5HNST AT . 9x I MS . 1 5X6HVEL0C Y/ ( I 6 * 2E 1 8 . 6 > ) 


CALL PRINT 0106 

GO TO 110 0107 

90 READ TAPE 8.J.MYIELD* C I NETO . ( R ( I ) . DR ( I ) . Z ( I ) . DZ ( I ) . DELTS ( I ) . OELTT 0108 

1 ( I ) *SINT( I ) .COST ( I ) . ( (SNT(L.K. 1 ) . SNP(L.K, 1 ) .Lnl ,NSFL ) .K= I .NFL) .1= 0109 

21 .Nl ) .TIME. (BIGNT ( 1 ) ,BIGNP ( I ) .B IGMT ( 1 ) .BI GMP ( I ) .YIELD ( I ) . I = 1 .N) 

CALL PRINT 
REWIND 8 

no BIGNTR(1)= 0. 0111 

QTR ( 1 ) = O . 0112 

DT(1)=0. 0113 

OTR(N+l )= 0. 0114 

DT (N+1 ) = 0. 01 15 

DS I NT ( N+ 1 ) = 0 . 0116 

DOT (N+1 )= 0. 0117 

120 J= J+1 0118 

D0S(1)3 2.*DR(1) 0119 

DDS(N+l )= (~2.)*DR(N) 0120 

CALL STRAIN 0121 

CALL STRESS 0122 

CALL EQUILI 0123 

OZ(N )«0. 

IF(J-M1 )I30. 140.210 0124 

130 IF(NYIELD) 150. 160. 150 0125 

140 Ml= M1+M2 0126 

150 TIME= DELTAT»FLOAT( J ) 0127 

CALL PRINT 0128 

NYIELD= 0 

160 IF (J-M ) 120. 1 70.210 0130 

170 1F(MYIELD)180.200.180 0131 

180 TIME= DELTAT*FLOAT (MYlELD) 0l32 

PRINT 190. time 0133 

190 FORMAT(28HO FIRST YIELDING AT TIME=E12.5> 0134 

200 WRITE TAPE8. J.MYIELD. C 1 NETO . ( R ( I > . DR ( I) . Z ( I ) . DZ ( I ) . DELTS ( I ) . DELTT 0135 
1 ( 1 ) .SINT ( 1 ) .COST ( 1 ) . ( (SNT(L.K. 1 ) * SNP (L.K. 1 ) .L = l .NSFL ) .K=l .NFD .1= 0136 

21 .Nl > .T IME. (BIGNT ( I ) ,BIGNP ( I ) .BIGMT( 1 > .BIGMP( I ) .YIELD ( I) . 1 = 1 .N) 

210 60 TO 5000 

END 0139 
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APPENDIX B - Continued 


Input Data 

The following data are input for a typical computer run for the clamped, 8-inch- 
diameter, circular plate. This particular run is for a = 0.5 inch, K2 = 1.5, 
gg = 1279 ft/sec, and h = 0.036 inch. The input to DEPROSS 3, as described in refer- 
ence 9, is given in the in-lb-sec system.; therefore, the input data presented in this 
appendix are given only in the in-lb-s system. 

The input data (see ref. 9) are as follows: 


N = 32 


NFL = 4 
NSFL = 2 


M = 18 000 
Ml = 200 
M2 = 200 
KE = 0 


= 1 

BIGR = 4.0 in. 

H = 0.036 in. 

CFL = 0. 

= 1 . 

SNU = 1/3 

RHO = 2.62012 X 10~4 Ib-sec2/in4 
DELTAT = 0.3 X 10~'7 sec 


number of mass points along the radius of 
sheet 

number of layers specifying the sheet 
thickness 

number of sublayers in strain-hardening 
model = number of coordinate pairs of 
stress and strain (i.e., 2) defining bilinear 
stress-strain curve 

cycle at which run is to stop 

cycle at which regular printing is to begin 

print every M2 cycles 

strains on the surfaces of the sheet are 
omitted in print sequence 

strains on the surfaces of the sheet are printed 

radius of circular sheet, b 

thickness of sheet 

layered cross section equivalent to fully 
plastic cross section 

layered cross section equivalent to fully 
elastic cross section 

Poisson's ratio 

mass density of sheet material 
time interval per cycle 


D, P = 0. pertain to constants used in strain-rate 

formula (see ref. 9) which were neglected in 
present application of DEPROSS 3 
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EPSIL (1) = 0.003627 
SIGMA (1) = 38080 lb/in2 
EPSIL (2) =0.15 
SIGMA (2) = 64064 lb/in2 
BIGA = 0.75 in. 

RHP = 0.80704143 X 10"4 Ib-sec2/in4 
GO = 15348 in./sec 
CIO = 1.5 

MV 

1 

0 


strain at yield stress 
stress at yield 

strain (approximate) at ultimate stress 

stress at ultimate strain 

radius of crushed hailstone, A 

mass density of hail (ice) 

hailstone velocity prior to impact 

multiplication factor related to the instanta- 
neous spread of the hailstone base at 
impact, K 2 

program instruction 

indicates an initial computer run and the 
velocity distribution is calculated by the 
program 

indicates a restart and the velocity distribu 
tion is read off a tape 


The input data appear on the data cards as 


3 ? 


1 


4 230000 

. 40000000+1 
. 2620 1 200-3 
. 36270000-2 
. 76 +0 


. 36000000-1 

.3 

.•^808 +5 

. 607041 43-4 


. 1 500+'000 + C 
. 1 6348000+5 


. 33333333+0 

.64064 +6 

. 1 60 OO 000+ 1 
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APPENDIX C 


CLASSICAL TfflN-PLATE BENDING ANALYSIS 


Determination of Radius of Central Plastic Region 

In reference 4 an elastic analysis was used to determine the radial and circumfer- 
ential bending-moment distributions in an infinite sheet impacted by a hailstone. The 
area of the sheet around the center of impact, where the moments exceed the elastic 
yield moment, was determined from the moment distributions of the elastic analysis. 
This central plastic region was then approximated by a simply supported plastic plate. 

A rigid-plastic analysis applied to this plastic plate determined the permanent deforma- 
tion or dent. 

The simple-support boundary conditions of the rigid-plastic analysis are Mj. = 0 
and Mq = Mq. Hence, the periphery of the central plastic region was taken to be at the 
radius | where Mj. = 0 and = Mq, simultaneously. This radius was determined 
in the following manner: By setting Mi.(t 7 ,t) equal to zero in equation (11) of 
reference 4 


DkVnl^ 


Mr(7?,T) ^ J3/2(p)sin(p2T)p"^/2 |3jj2(p7j) - JoCprjJ - | p-lji(p7?)j dp = 0 (Cl) 


specific values of the nondimensional radial location rj can be determined for various 
values of the nondimensional time parameter t. With these specific values of tj = ^ 
and T and equated to Mq in equation (12) of reference 4 


DkVnl^ 


Mg = 


6A^ 


J J3/2(p)sin(p2T)p"V2 


■f + J2(P^) 


Jo(p^) 


dp = Mq 


(C2) 


corresponding values of hailstone velocity ^Vq or ggj can be determined for particular 
sheet and hailstone properties of interest. 

In figures 11 and 12, respectively, the dimensionless radius of the plastic region | 
is plotted as a function of hailstone velocity gQ for a 0.122-cm- (0.048-inch-) thick L 72 
and L73 aluminum sheet. A curve is given in each figure for the 1.27-, 1.91-, and 
2.54-cm- (0. 5-, 0. 75-, and 1-inch-) diameter hailstones (represented in the figures by 
their radii). In the calculation of these curves, the hailstone radii were increased by the 
appropriate K 2 factor to account for the experimentally observed spreading discussed 
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APPENDIX C - Concluded 

previously. The radius of the plastic region and other pertinent sheet and hailstone 
parameters can now be used to determine the permanent deformation or dent. 


Determination of Permanent Center Deformation 

The permanent center deflection 6 of the plastically deformed plate can be deter- 
mined from the following equation (see eq. (16) of ref. 4): 


6 = 


MV§a2|2 


12Mf 


1 4 \ 




r 2^ 


_3 

H 


3 "'"“' I 


(1 S 4 S oo) 


6 = 


MV§a2^2 


12Mr 


.1+-L+ Jl - |2/i _ J 
2 ^2 'I I 2^2 


2|‘ 


sin- 1(1) 


(I S 1) 


(C3) 


with the appropriate values of A, p, h, Mq, Pg, Pp, g^, and Vq, where 


Vn = 


h Ps 
h Ps 


(C4) 


and il/ = 0.70. 
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TABLE I.- MATERIAL PROPERTIES^ 


Material 

Yield stress 

Yield strain 

Ultimate stress 

Ultimate 

strain 

— 

Density 

Modulus of 
elasticity 


psi 

MN/m2 


psi 

MN/m2 

Ib-s2/in4 

kg/m3 

psi 

GN/m2 

L72 (2014-T4) 

38 000 

262 

3.63 X 10-3 

64 000 

441.3 

0.150 

2.62 X 10-4 

2800 

10.5 X 106 

72.4 

L73 (2014-T6) 

58 250 

401.6 

5.55 X 10-3 

67 000 

462 

.090 

2.62 X 10-4 

2800 

10.5 X 10® 

72.4 

Ice 







.81 X 10-4 

866 




^The input to DEPROSS 3, as described in reference 9, is given in the in-lb-s system. 


ts:) 

w 





(b) Velocity and/or mass distribution approximations used in DEPROSS. 

♦ 



(c) Initial velocity distribution of clamped plate and crushed hailstone. 


Figure 1.- Mass distribution, velocity and/or mass distribution approx- 
imations, and initial velocity distribution for a clamped circular 
plate impacted by a hailstone. 
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Scale; 


1 inch (2. 54 cm) 



(a) Normal impact; a = 0.635 cm (0.25 inch); go = 658 m/s (2160 ft/s); 

h = 0.091 cm (0.036 inch). 



(b) Normal impact; a = 0.952 cm (0.375 inch); go = 547.1 m/s (1795 ft/s); 

h = 0.122 cm (0.048 inch). 

Figure 2.- Computed shapes at different times for two specific hail impacts. 

(Deformations to scale.) 
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Time, us 


Figure 3.- Center deflection as a function of time for a 1.27-cm- (0.50-inch-) diameter hailstone 
impacting a 0.091-cm- (0.036-inch-) thick L72 aluminum sheet at various velocities. 


rrnrr 


Scale: 


inch (2.54 cm) 

- DEPROSS 

- R.A.E. date 



(a) Normal impact; a = 0.635 cm (0.25 inch); h = 0.091 cm (0.036 inch). 



(b) Normal impact; a = 1.27 cm (0.50 inch); h = 0.091 cm (0.036 inch). 

Figure 4.- Comparison of theoretical and experimental sheet deflections for 
L72 aluminum. (Deformations to scale.) 
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Figure 6.- Comparison of theoretical 
impacted by 1.91 
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Figure 12.- Radius of plastic region as a function of hailstone velocity for aluminum 

alloy L73. h = 0.122 cm (0.048 inch). 
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